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The case of gluon bremsstrahlung off a heavy quark in extended nuclear matter is revisited within 
the higher twist formalism. In particular, the in-medium modification of “semi-hard” heavy quarks 
is studied, where the momentum of the heavy quark is larger but comparable to the mass of the 
heavy quark (p > M). In contrast to all prior calculations, where the gluon emission spectrum is 
entirely controlled by the transverse momentum diffusion parameter (g), both for light and heavy 
quarks, in this work, we demonstrate that the gluon emission spectrum for a heavy quark (unlike 
that for flavors) is also sensitive to e, which so far has been used to quantify the amount of light-cone 
drag experienced by a parton. This mass dependent effect, due to the non-light-like momentum of a 
semi-hard heavy-quark, leads to an additional energy loss term for heavy-quarks, while resulting in 
a negligible modification of light flavor (and high energy heavy flavor) loss. This result can be used 
to estimate the value of this sub-leading non-perturbative jet transport parameter (e) from heavy 
flavor suppression in ultra-relativistic heavy-ion collisions. 

PACS numbers: 12.38.Mh, 12.38.-t, 12.38.Cy, 12.38.Bx 


I. INTRODUCTION 

The unprecedented centre of mass energies available 
at the Large Hadron Collider have opened new windows 
for the exploration of extreme nuclear matter through 
high energy jets While a large portion of the 

available data on leading (and next-to-leading) particle 
suppression in the light flavor sector has been theoreti¬ 
cally described using factorized pQCD based calculations 
of jet modification [3, 0, heavy quarks have remained 
somewhat of a challenge [^. This is especially true in 
the semi-hard sector, where the momentum of the heavy 
quark is larger but comparable to its mass p > mg. We 
distinguish this region from that of slow heavy quarks, 
where p ^ mg, which appear to be thermalized with 
the bulk medium, and fast heavy-quarks, with p ^ mg 
which engender energy loss and suppression similar to 
light quarks. 

The so called “heavy-quark puzzle” had already be¬ 
gun to manifest itself in measurements of the suppression 
of high transverse momentum (high-px) non-photonic 
electrons at the Relativistic Heavy-Ion Collider (RHIC). 
Measurements by both the STAR [l^ and PHENIX [h| 
detectors showed a slightly higher suppression than ex¬ 
pected, based on a calculation that included both drag 
and radiative loss [1, [12, [12] • This trend has continued 
at the Large Hadron Collider (LHC) where the ALICE 
experiment has measured D and B meson suppression 
separately, and finds a larger than expected suppression 
in the semi-hard regime of heavy-quark momentum (we 
note that the case is not very clear for B-meson suppres¬ 
sion which has so far only been presented as pt integrated 
points) [II m. 

A considerable amount of theoretical work has been 
devoted to understand this larger than expected suppres¬ 
sion of single electrons or heavy mesons arising from the 
fragmentation of a heavy-quark. However, most of these 
may be understood as falling in two categories: Calcula¬ 


tions that have extended the base formalism of radiated 
energy loss for light flavors to include mass dependent 
terms, as well as a drag term to include the prominent 
role played by drag in heavy flavor energy loss [i,[il- 
[^ . Calculations that have ignored the role of radiative 
loss and only focussed on drag loss [ 2 ll - [^ . 

In all calculations above, radiative loss is stimulated 
by transverse momentum diffusion experienced by the 
heavy quark or radiated gluon, which, in some cases, is 
quantified by the jet transport coefficient q [H, [2^. The 
drag loss is quantified using the drag coefficient referred 
to as dE/dx (energy loss per unit distance) or e [ 1 ^ . 

To the best of our knowledge, no calculation of heavy 
flavor energy loss has explored the possibility that the 
drag coefficient e (or the longitudinal diffusion coeffi¬ 
cient 62 ) rnay lead to an additional source of radiative 
loss, beyond that provided by q. This possibility is im¬ 
mediately clear in the higher twist framework, where the 
drag (and longitudinal diffusion) coefficient e ( 62 ) has the 
boost invariant definition as the loss of light-cone mo¬ 
mentum (fluctuation in light-cone momentum) per unit 
light-cone length, (assuming a parton moving in the neg¬ 
ative light-cone direction) 

. d{Ap-} . d{Ap-^) 

dL- ’ dL- ' ^ ^ 

While such a transport coefficient leads to little change 
in the off-shellness of a near on-shell massless quark, it 
has a considerable impact on the off-shellness of a near 
on-shell massive quark. No doubt, such a term will only 
have an effect on the radiative loss of a patron where the 
momentum p is comparable to the mass M, thus for light 
flavors, and for energetic heavy-quarks where, p ^ M, 
this term will have a minimal effect. This was explic¬ 
itly explored for photon radiation from a light quark in 
Ref. [ 23 ]. We point out that such an effect is by no means 
limited to the higher-twist scheme, but effects several 
other formalisms that have considered the radiative loss 
from a heavy-quark in a quark gluon plasma. 
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In this paper, we will explore the modification to the 
calculation of radiative loss, due to the presence of this 
additional source within the higher twist formalism. To 
delineate the importance of these terms, we will use 
power-counting techniques borrowed from Soft-Collinear- 
Effective-Theory (SCET) [ 23 - 1 ^ to identify the regime 
where these mass dependent terms will cause detectable 
effects on the gluon bremsstrahlung spectrum. As a first 
attempt, we will consider only the case of single scat¬ 
tering and single emission. Arguments presented in sub¬ 
sequent sections will demonstrate this to be the leading 
contribution, given the short formation time of the ra¬ 
diated gluons. In this paper, the analytical expressions 
for the e {q and 62) induced gluon radiation spectrum 
will be derived. Numerical calculations for the suppres¬ 
sion of B and D mesons, as well as the suppression and 
azimuthal anisotropy of non-photonic leptons from the 
decay of these mesons at LHC and RHIC energies, and 
comparisons with experimental data will be carried out 
in a subsequent effort. 

The article is organized as follows: in Sec. II, we will 
setup the basic formalism of Deep Inelastic Scattering 
(DIS) on a large nucleus, where the hard virtual photon 
strikes a heavy-quark, assumed to be produced in a rare 
high fluctuation inside a proton. In Sec. Ill, we will 
carry out diagrammatic studies on the induced gluon ra¬ 
diation off the heavy quark in this system, within the 
higher twist formalism. We will present and discuss final 
expressions for gluon bremsstrahlung from such a heavy 
quark, containing both g, e and 62 in Sec. IV. We offer 
concluding discussions and an outlook in Sec. V. Involved 
expressions for a set of diagrams are contained in the ap¬ 
pendix. 


II. DEEP INELASTIC SCATTERING AND THE 
SEMI-HARD HEAVY QUARK 

The set-up is based on the deep-inelastic scattering of 
a virtual photon off a heavy quark within a large nucleus 
with mass number A. We will study the case where the 
hard virtual photon scatters with the hard heavy quark 
converting it to slow moving heavy quark (this is defined 
below). The propagation of the heavy-quark will be fac¬ 
torized from the hard scattering vertex which produces 
the outgoing slow moving heavy-quark. The nucleus has 
a momentum P = pA, where p is the average momentum 
of a nucleon in this nucleus. In the Breit frame, the ex¬ 
changed virtual photon possesses no transverse momen¬ 
tum, 

q=[q~'',q~,qi-]= [q^,q~,o]. ( 2 ) 

The process under consideration is the following: 

b{Li) -f A{P) —>■ e{L 2 ) -f Jq{Lq) X. ( 3 ) 

In the above process, e(Li) and 6(^2) represent the in¬ 
coming (outgoing) electron with momentum Li and L2 
respectively. The factor A{P) represents the incoming 


nucleus with momentum P. The factor Jq{Lq) is the 
outgoing jet which contains one heavy quark Q. In 
the particular diagrams that we will consider, the final 
strongly interacting particles Due to the absence of va¬ 
lence heavy-quarks within the nucleon, the photon will 
have to strike a heavy quark from s QQ fluctuation within 
the sea of partons. Alternatively one may consider the 
case of a, J/tjj ov T state bound within a large nucleus, 
being struck by the hard virtual photon. More detail dis¬ 
cussions on the production of heavy quark can be found 
in [s^. In this work, we will not discuss the produc¬ 
tion of the heavy quark further. For the purposes of this 
calculation, this is now contained within a parton distri¬ 
bution function. In essence, by this mechanism a semi- 
hard heavy quark has been produced. In what follows, 
we will focus on the power counting of the momentum 
components and the modification of the final state. 

Similar to Ref (3^ . we will consider a quark mass 
M ^ Aqc u and a final outgoing quark momentum which 
is larger, but of the order of the quark mass. We assume 
that the quark, anti-quark fluctuations possess minimal 
momentum in the rest frame of the nucleus, and thus 
its momentum components scales as (pq,Pq,Pq_l) 
(M/-\/ 2 , M/^/ 2 , 0). Now in a frame where the nucleus 
is boosted by a large boost factor 7 in the “-I-” direction, 
momentum components of the incoming heavy quark will 
scales as. 


PC = [pc>Pq:PC-l] = 


' M I M ^ 

771 ’ 


( 4 ) 


It is important to note that the boost factor 7 carries 
no extra information other than the the fact that p'q is 
very large compare to pg. We assume the momentum 
components of the incoming photon to be, 


q = 




I M 

771 ’ 


( 5 ) 


Given a large boost factor (7), one can assume that 
7M ^ M ~ ^ M/7. Hence, we define the hard scale 

Q as = —q^ ~ ^Mq~/\/ 2 . As a result, we obtain the 
final out-going quark to have momentum components as. 


Pf =PQ + q 


M2 

2 ^ 


,q ,0 


( 6 ) 


We consider M < q~ ^ \f\Q where A ^ I/7 for this 
“semi-hard” heavy quark. The term “semi-hard” defines 
a quark whose momentum is not an order of magnitude 
larger than its mass. 


A. Power Counting and the small A parameter 

In order to set up the power counting, in this study, 
we have introduced the dimensionless small parameter 
A. Power corrections to hard process are generally sup¬ 
pressed by factors of a hard scale, ^ Aqcd- The 
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introduction of the parameter A to represent semi-hard 
scales as XQ and softer scales as X^Q, is a concept bor¬ 
rowed from soft collinear effective theory (SCET) [s^, 
Is^l - In what follows, we will retain leading and next-to- 
leading terms in A power counting, neglec ting all terms 
which scale with A^ or a higher power of A . We have 
chosen the scaling variable A in such a way that perturba¬ 
tion theory may be applied down to momentum transfer 
scales at or above A^/^Q ~ Aqcd- In our study, 

A°>VA»A»A5. (7) 


Based on the power counting set up and the choice of 
incoming and outgoing quark momentum, we can outline 
several important scales in the problem of a heavy quark 
propagating through the nuclear medium, scattering off 
constituents in the medium and emitting real gluons. We 
remind the reader that given the semi-hard momentum 
of the heavy quark, collinear emission is suppressed due 
to the mass of the heavy quark, i.e., the dead cone ef¬ 
fect [ 3 ^. In the subsequent section, the calculation of 
the process of single scattering and single emission from 
a heavy quark will be carried out. Here we outline the 
power counting (in A) of the relevant momentum compo¬ 
nents that will arise in the calculation. The virtuality of 
the hard photon defines the hardest scale in the problem, 
(5, similar to the case of light quark production in DIS. 
The incoming or initial heavy quark has momentum com¬ 
ponents Pi ^ (A “2 , A 2 ,0)(5, the outgoing heavy quark 
has momentum components p ^ (V^, a/A, 0)Q- The mass 
of the semi-hard heavy quark scales as M ~ a/AQ- This 
choice of incoming parton and photon momenta ensures 
that the momentum of the final outgoing heavy quark 
is of the order of its momentum. In what follows, we 
will demonstrate that the leading contribution to gluon 
emission arises from the region where real emitted gluons 
have momenta which scale as I ~ (A, A, X)Q. The frac¬ 
tion of light cone momenta carried out by the gluon is 
y = l~/p~ ~ VX. Also somewhat different from the case 
of light flavors is the scaling of the virtual Glauber gluons 
: k ^ {X^,Xi,X)Q with = 2k'^k~ — ~ As 

we will demonstrate below, these choices of momentum 
scales tend to enhance the gluon emission rate. 

There is another consequence of this choice of scales 
which relates to how the single gluon emission kernel 
may be iterated. Unlike the case for light flavors, the 
formation length of a gluon with momentum components 
I ~ (A, A, X)Q is 



1 

XQ’ 


( 8 ) 


is rather short compared to the formation length of 
Tq ^ 1/X^Q for gluon radiation from near on-shell light 
flavors. This indicates that there cannot be many scat¬ 
terings per emission. As a result, in what follows, we 
derive the single scattering per gluon emission rate. This 
single gluon emission kernel, induced by single scattering 
will have to be iterated to obtain the full energy loss of 
a semi-hard quark. 


Many readers may find the presence of factors of y/X 
somewhat disconcerting. We could have simply replaced 
this with a new A. We refrain from defining a new dimen¬ 
sionless parameter A, so as to make contact with prior 
dehnitions of A used in the case of light quarks, where 
XQ represented the transverse momentum of the radiated 
gluons from a hard parton, or the transverse momentum 
from scattering of a gluon in the medium. Thus, to con¬ 
tinue to draw a parallel with the prior results form light 
flavor energy loss, we require XQ ^ 1 GeV. 

To get a physical feel of these scaling relations, one may 
typically assume Q ~ 100 GeV, y/XQ ~ 10 GeV, XQ ~ 
1 GeV, X^Q Aqcd- We will retain terms that are 
0{-\/X) suppressed compared to the leading terms, but 
ignore all terms that are suppressed by 0(A) and higher. 
Thus, terms with ~ A will eventually be ignored. 


III. INDUCED GLUON RADIATION OFF THE 
HEAVY QUARK 

In this section we will discuss some contributions to the 
next-to-leading order correction to semi-inclusive DIS on 
a large nucleus with a quark and gluon in the final state. 
By next-to-leading order, we simply mean including one 
interaction term in the amplitude and complex conjugate, 
which converts a single quark to a quark and a gluon. 
The diagrams we will consider will also contain two scat¬ 
terings for the full cross section, with both scatterings 
in the amplitude, both in the complex conjugate, or one 
in amplitude and one in complex conjugate. The double 
differential cross section for the semi inclusive process of 
an electron with in-coming momentum Li and out-going 
momentum of L 2 , undergoing DIS off a nucleus (with mo¬ 
mentum pA), leading to the production of a final state 
heavy quark with transverse momentum Iq^ and a final 
state gluon with transverse momentum ^ Xlqcd, may 
be expressed as 


Eh^da ^ j dW^'' 

d^L2d?'lQ^d?l±dy 2 ttsQ'^ d^lQ^d^l^dy 

In the equation above, s = (p-\-Li)^ is the total invariant 
mass of the lepton nucleon system. In the single photon 
exchange approximation, the leptonic part of the cross 
section is easily expressed in terms of the leptonic tensor 
denoted as given as, 

Dpi/ = iTr[,£i7^ ^2lu\- (10) 

In what follows, the focus will lie entirely on the 
hadronic tensor W^'^. We will carryout calculations 
of a set of contributions to at next-to-leading or¬ 

der (NLO) as described above, and next-to-leading twist 
(NET), meaning double scattering in the cross-section. 

Already, at NLO and NLT, there are several interfer¬ 
ing diagrams to consider. In this section, the calculation 
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FIG. 1. A representative single gluon emission diagram, 
where gluon emission is induced by single scattering.This rep¬ 
resents a symmetric diagram with scattering off the final pro¬ 
duced quark. Three separate cuts, denoted as central, left 
and right are indicated by the dashed lines. 


of one of the diagrams that contribute to single scatter¬ 
ing induced single gluon emission will be carried out in 
some detail to familiarize the reader to the approxima¬ 


tions carried out in this article. Figure o represents 
the diagram that will be evaluated. This diagram cor¬ 
responds to the process where a semi-hard heavy quark 
produced after DIS radiates a gluon followed by a scat¬ 
tering in the medium, and finally exits the nucleus. 

In this study, calculations will be carried out in axial 
gauge, n ■ A = 0, with n = (1,0, Oj.) and A~ = 0. In 
A~ = 0 gauge, the double scattering of a quark radiat¬ 
ing a gluon contains a total of nine central cut diagrams 
where the cut line lies between the two scatterings. It also 
contains seven left cut diagrams and seven right cut dia¬ 
grams. In this section, the central cut diagram of Fig.[T] 
will be analyzed in detail. In this section and what fol¬ 
lows, only the real contribution where the radiated gluon 
line has been cut will be considered. The entire contribu¬ 
tion from virtual diagrams, which contain a quark gluon 
fluctuation either in the amplitude or complex conjugate, 
will be obtained using unitarity arguments. 

This NLO-NLT contribution to the hadronic tensor 
may be expressed as. 


= 9\-gT) 


■ d%d^p', d^d'^k' d^yid^k d^y^d^p, dH dHg ^-ipivo 


(27r)4 (27r)4 (27r)4 (27r)4 (27r)4 (27r)4 


xTr 


i\4-2 ;^2 A/r2 ;^2 H/r2 | I\2 1\4^2 I 


2 "^ p'j — — ie {pf — ly — {pf — lY — -|- ie {pf — lY — + ie 


x2TTS{Y)27rS{YQ - M^){-iT^){tT^){A\Y{y'Yl^A+{y[)A+{yYY{yo)\A). 


( 11 ) 


We have defined the following momentum fractions for 
convenience, 


I k l-y 

q I l-y+ r]y 

pf k+ 


X = 


XL = 


XK = 


‘j. 


2P+q-y{l-y) ’ 


xd = 


kl 


2P+q- ’ 


K = 


1 + (1 - y) 


y^M^ 

~P~ ’ 

k\ — 2lj_k± 
2P+q- 

r , xm = 


( 12 ) 


IVP 


2P+q- 


In what follows, we will simply both numerator and de¬ 
nominator by retaining only leading terms in the A power 
counting highlighted above. This is followed by taking 


the trace in the numerator and contour integrations to 
simplify the denominator. 

The factor is the polarization tensor, and in 

A~ = 0 gauge, it has the form, 


Qay _ _gO‘P 


nl 


(13) 


The structure of n = (1,0, 0 _l) implies n-l = l~. Within 
this kinematic set up, in A~ = 0 gauge, the leading com¬ 
ponent of the gluon field is A+. 

The spin sum in the numerator, containing the entire 
set of 7 matrices from the quark propagators, together 
with those from the interaction with the gauge field may 
be partially simplified and expressed as. 




























5 


ATf 1 = Tr 


27 7 9 7.- 


^ii±y 


■7+ (g -I ) + 7“p/±„7-Lm (-ffD 7"^ (g - O 


+7’^g 7-LM(-gr)7“ (p/i-^- l)„} 7 7’^(g -i +fc )7 |7’^(g -i )ip —p 

+7’^(g" - ^■)7l(-g±)p^7±p/_L;3 + T'± {Pf± ~ ^i);3 7l(-g-L)p!/7’^g”}] 

'1 


+ _ 

-7 g 


+ Tr 


-1 {M 7 _l^ (-gr) 7“^ (g - O + 7-Lp (-^r) 


X 7 7 ^(g +fc )7 { 7 +(g -I ) 7 ^(-g_L)pj,M + M 7 ^(-g_L)pi. 7+9 }] . 


(14) 


Note that terms containing 7 pj never contribute to the 
trace because there is always an adjacent factor of 7“, 
and 7“7“ = 7+7+ = 0 . One may evaluate the trace, 
using the relation 7“''7“ = 2 — 7“7^ Ea. dHl) simplifies 
to, 


J^c ^ ‘^{‘^<1 ) (1 _ y + yj^) [P{y)l\ + . (15) 

y 


Note that the mass independent portion contains the 
standard vacuum splitting function P{y) while the mass 
dependent part has a separate dependence on the mo¬ 
mentum fraction of the radiated gluon (y). In the soft 
emission kinematic limit where y <C 1, one may neglect 
the mass term. However, in this work we will retain it 
throughout. 

The set of denominators can now be evaluated us¬ 
ing contour integration. For the central cut diagram of 
Fig.O, this yields the phase factor, 


III = exp 


l \ Xb + XL + 


Xm 


1-y 

X 6 '(yr -Po") ^ [y'l -g'o) 


P^ (y' 


0 ~ Po 


exp 


• ( A , /A A W 

* Qxd -f (C - 1 ) -7 - C:;- -XL 


1 - y 1 - y 


P^ (y'l 


yi 


1 — exp 


-i[xL + “^0) 


1 — exp 


i\XL + 


1-y 


XM P 


(y'l -y'o) 


This single diagram contains four contributions depend¬ 
ing on the one propagator that remains off shell after 
contour integration. 

Multiplying through we find four terms similar to the 
work of Ref. [ 37 j . The first term where the same propaga¬ 
tor is off-shell in both amplitude and complex conjugate 
corresponds to the so called “hard-soft” process where 
the gluon radiation is induced by the initial hard scatter¬ 
ing. The heavy quark is knocked off-shell by the initial 
hard scattering and becomes on-shell after radiating the 
on-shell gluon. Afterwards, the on-shell quark or gluon 
will have a scattering with another soft medium gluon 
from the nucleus. The second term is the case where the 
quark is on-shell immediately after the first hard scatter¬ 
ing. Gluon radiation is induced by subsequent scattering 
of the heavy quark off a in-medium gluon which carries a 
specific finite momentum fraction. This is often referred 
to as “hard-hard” scattering. The two cross terms where 
different propagators are off-shell in the amplitude and 
complex conjugate represent interference between soft- 
hard and hard-hard scatterings. 

The equations derived above contain both longitudinal 
and transverse momentum exchanges with the medium. 


The portion due to transverse exchange may be isolated 
by imposing that —>■ 0 (y —>■ 0, C 1) limit, and then 
comparing with expressions from similar diagrams in Ref 
(2^. One will immediately note that factors contain¬ 
ing Xm which contribute to the Landau-Pomeranchuck- 
Migdal (LPM) [H, [s^ effect are not modified by factors 
of while factors containing xd are modified by the 
presence of Factors oixu will eventually be absorbed 
in the definition of the transport coefficients including q. 
Such factors introduce a non-trivial dependence of in¬ 
medium transport coefficients on the mass of the probe. 


In this section, we have demonstrated how a certain 
diagram for heavy quark production and energy loss via 
gluon radiation can be simplified. Similar rules will be 
applied to all other real diagrams which include a cut of 
the radiated gluon line. In the subsequent section, all 
real diagrams will be combined to obtain the real gluon 
emission spectrum from a heavy quark that undergoes 
one scattering and one emission after production. 
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IV. GLUON EMISSION SPECTRUM 


In the preceding section we evaluated the diagram in 
Fig. [ 1 ] in some detail, to highlight the approximations 
that will be made in the course of the full calculation. 
In this section, the result of the sum of all real diagrams 
(with an emitted gluon in the final state), will be pre¬ 
sented. This will be followed by a gradient expansion in 
the exchanged transverse momentum (fcj_ —> 0 ). While 
the leading term in the limit of fcj_ —> 0, will correspond 
to a gauge correction to the vacuum process of gluon 
radiation from a heavy quark, the focus in this section 
will be on first correction in the A:_l —>■ 0 limit, usually 
denoted as the next-to-leading twist contribution. 

In total, there are II separate topologically distinct di¬ 
agrams similar to that in Fig. [T] We denote these with 
two subscripts: m,n = 1 - 3 , where, m denotes the loca¬ 
tion of the scattering in the amplitude, and n denotes 
the location in the complex conjugate, for the case of a 
central cut, where one gluon scattering is on either side 
of the cut line. In either case of to or n, 1 signifies that 
the scattering occurs on the quark line beyond emission, 
2 signifies scattering on the quark line between the hard 
production and the emission, and 3 signifies scattering of 
the emitted gluon. Each one of these diagrams will also 
generate a left and right cut component, where the cut 
line will be moved to the left or right of the scatterings, 
with the topology of the diagram held fixed. There are 
also the two special configurations, where both scatter¬ 
ings occur between the hard production and the gluon 
emission in the amplitude or complex-conjugate. These 
are denoted as Co,i and Ci^. The next-to-leading twist 
portion, of the sum over all three cuts, for each of these 
contributions is outlined in the appendix. 

Adding all the contributions from all the diagrams, 
categorized in the Appendix, we obtain the entire con¬ 
tribution to the hadronic tensor. In what follows, we 
decompose the hadronic tensor as. 


( 16 ) 


The entire contribution from all real diagrams is con¬ 
tained in the factor This, includes the initial hard 

scattering, the final state scattering of the quark or gluon, 
and the emission vertex. Virtual contributions, where the 
final state radiated gluon is not cut, will not be consid¬ 
ered in this effort. Some part of the spin sum in the 
numerator has already been factorized out in the term 
—g^ above. In what follows, we will simplify by 

factoring different contributions within it and then apply¬ 
ing approximations to them separately. In the interest of 
readability, the exact details of the calculation for each 
diagram separately is included in the appendix. 


This entire factor is obtained as, 

3 

_ \ ' nc,l,r,ri , n _ 

n — 2_^ -|-Co,i-l-Ci,o — 

m,n=l 


(y'o - Vo ) 

+ *(cx,, + (C-i)^-C^Xi) 

V 1 -y 1 -y J 

X P"" {y'l -yi) 

X {A\ij{y'Q)-i+A+{y[)A+{yi)i){yo)\A) (17) 

In the equation above, cut specific phase factors and 
the hard part for each cut are entirely contained within 
overall phase factor represents the generic 
portion of the phase factor. 

In order to calculate the next-to-leading power contri¬ 
bution to the semi-inclusive hard partonic cross-section, 
one needs to expand the cross-section in fe_L and in k~. 
In each case, we will extract the corresponding trans¬ 
port coefficients inside the gluon emission spectrum ker¬ 
nel for the semi-hard heavy quark . Factors of fcj_ and 
k~ are absorbed as derivatives within the definition of 
the transport coefficients [e.g. k±A'^{y±) exp(ffcj_ • y±) = 
-iV±A+ {y±) exp{ik± ■ y±) ~ iF+-'-(y±) exp(ikj_ ■ y±)]. 

We will also factor the four point non-perturbative op¬ 
erator using the usual phenomenological factorization, 
which for the case of transverse scattering may be ex¬ 
pressed as, 

{A\^{yoh^F+{y[)F+{yi)^{yo)\A) ~ 

X (p|V'(yo)7^V'(2/o)b)^(p|A+(yi)F+(yi)|p). (18) 

The first operator product on the right hand side of the 
equation above will yield the incoming quark distribu¬ 
tion function within one nucleon. The second operator 
product will yield the transport coefficient due to the 
scattering of the final state, off a gluon within a nucleon 
in the nucleus. We have assumed the average condition 
that both nucleons have a momentum p = P/A. The fac¬ 
tor p represents the nucleon density within the nucleus, 
and Cp represents an overall normalization constant that 
contains the nucleon density. The factor of p/( 2 p+) is 
written separately as that will be absorbed within the 
definition of the transport coefficient. 

These transport coefficients, defined below, are non- 
perturbative objects, which are factorized from the hard 
part that describes the propagation of the heavy quark. 
While the exact value of each of these coefficients de¬ 
pends on the non-perturbative dynamics of the medium, 
the relative contributions of the different hard parts that 
appear as a multiplicative factor along with these coeffi¬ 
cients will be calculated below. Terms for the transverse 
diffusion q, the drag (and longitudinal diffusion) coeffi¬ 
cient e (62) can be obtained through derivatives of the 
kernel with respect to the transverse and (—)-light-cone 


X exp 


XB+XL + 


Xm 

1 — 7 / 


27 rQ:s 
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component of the exchange momentum, 
'l + (l-y) 2 ' 


Vfe-, Vfc- 


I fej_ ,fe“ =0 


= ^Ca 


y 


‘L 


[^i + y^M 


nr, 


214 C,i,r 


(19) 


In the equation above, the factor represents 

several terms, depending on the cut taken i.e., central 
c, left /, or right r, and the momentum component with 

I 


respect to which the Taylor expansion is carried out, i.e., 
q for the second derivative in terms of fcj_, e for the 
first derivative with respect to k~ and 62 for the sec¬ 
ond derivative with respect to k~. One should note that 
for each case, once the derivatives have been taken, both 
factors of the momentum k±,k~ are set to zero. The 
complete expressions for can be expressed as a 

sum of products of a phase factor and a non-phase factor 
coefficient, expressed as 


^ _ ^-i{xL + J^XM)P*{vi -Vg ) I _ + -y'o ) _|_ ^^-i{x:L + J^XM)P'*'{yi -yg ) 

X I - _|_ I _ Q-iixL+Y^XM)P+{yi-yo) ^i{xL+T^XM)P'*'{y'T-y'g)'^ 








^-i[xL + j^xM)P^{yZ-yo ) gi{xL + T&^XM)P^{y'i -y'o ) 

^-i(^XL+Y^XM)P+{y'o-y'i) _ g-i{xL+j^xM)P+{y'o-yT) 
-i(^xL + j^xM)P^{yo-yZ) _ p-i{xL + T^xM)P^{yo-y'i ) 




1 _ g-*(®i + T^a:M)P+(y'o -y'l) 




1 _ g-■i■{^L+l^xM)P+{yo -yi ) 


In each case above, the coefficients depend on the 

momentum component being considered. The subscript 
n merely denotes the order in which the coefficient occurs: 
Cl, C 2 and C 3 appear in the expression for the central 
cut, where C 4 and C 5 appear in both the left and right 
cuts. We list them in the following for each different 
case, starting from the case of transverse diffusion, i.e., 
q. The coefficients are. 


cf = 1 - (2 - 3Ky^)x + (1 - Ky )x , 


q y 

^ 2 


1 - -ny^) + y^^ (2 - ny^) 


X 


^ {y-P'y^)-y^^ i^y^) 


4 = [1 - 4(1 - Ky4x + (1 - 2«;y^)x^] 


^9 — 


1-f 


1-f 


y 


i + (i-y)^ 

,2 


y 


l + (l-y)2 


X 


X 


+ 1 - 2y 


X ■ 


( 20 ) 


The momentum fractions y, k and x ^^e defined in 
Eq. (fT^ . For the longitudinal drag coefficient e, the c- 

I 


factors are, 


y^M^ 

l- 

y^h/P 

l- 


- + - {1 + Ky 4 x+ 


1 1 I'l I 2\ 2 1 2 

- ^(1 + Ky ) y X - ^x 


c^ = 0 ,c^ = 0 ,c® = 0 . 

For longitudinal diffusion coefficient 62 , the c-factors are. 


y'^M'^ 

r 1 


1 2 

^ , 722 ] 

[l-f 

."2^ 

(2 


\x + i^y X 

y'^M^ 


3 

1 2\ 

3 2 2 

{l-f 

[4 ( 

4 

ry) 

x-^p-y X 

0 , c®^ 

= 0 , 

c®^ 

= 0 . 



All the terms presented above can be combined to ob¬ 
tain the real single gluon emission spectrum. In the sec¬ 
ond line of the equation below [Eq. (I?T]) ]. we have re¬ 
tained terms only up to the approximation that 

has been justified in this study. All terms which scale 
as 0{X) or greater, have been neglected. We express the 
gluon spectrum per unit light-cone length as. 


= 2 - P(v) - _ 

dydl^dr tt (^l-l-x^ 

= 2 ^P{y)^f^ 
^1 Vl + X 


(t^) ( 4 /-(i^-y) (^+x) [{4 + 4 } g + 4{cf + c®} e + 2{c?+c®^} 


62 


Sin 


:(1 +x) 


4/-(l - y) 

{(i-|)-x+(i-|)x"} «+bx(i + x)’ ' + 7pyx(I-)Ix)fe 


X 


( 21 ) 





























































In the equation above, we have defined a mean light-cone 
location of the first scattering (between the amplitude 
and complex conjugate) as 


Vi 


y'l 


( 22 ) 


We also define the off-set between the light cone locations 
in the amplitude and complex-conjugate as, 


y =yi -y’l- (23) 

This variable enters the definitions of all transport co¬ 
efficients that will be discussed in this paper. There 
are three transport coefficients, which contain the non- 
perturbative expection of the gluon field strength oper¬ 
ators: the transverse diffusion coefficient, q which repre¬ 
sents the transverse momentum squared per unit light- 
cone length, exchanged between the hard quark and the 
medium, the longitudinal drag per unit light-cone length, 
e caused due to the exchange of light-cone components 
of momentum, and 62 the diffusion in light-cone momen¬ 
tum, per unit light-cone length : 



An'^Cnas 

f dy 

P 


9 — 

iV2-l j 

1 7T 

2p+ 


X 

{A\F+{y- 

)F^+{0)\A) e- 



A'K'^CRas 

f dy 

P 



Nl-l J 

1 TT 

2p+ 


X 

{A\id~{y~)A+ {Q)\A) e 

—iKP~^y~ 



f dy 

P 


62 = 

iV2-l j 


2p+ 


X 

{A\F-+{y- 

-)F-+ 

(0)|A) e- 

-zAP’*'y~ 


In the equations above, we observe the appearance of 
another momentum fraction: 

A = Cxz. + (C - I) ^ . (25) 

1 - y 1 - y 

The presence of such a momentum fraction, indicates 
that the range of momentum fractions in the definition 
of q e and 62 for heavy quark scattering is different from 
that for light flavor energy loss. This indicates that the 
actual value of q (or even e or 62 ) for heavy quarks may 
be different from that for light quarks. Thus, careful 
analysis of heavy-quark energy loss may lead to an un¬ 
derstanding of the cc-dependence of the in-medium gluon 
distribution function that sources transport coefficients, 
and may, in the end, lead to an understanding of the de¬ 
grees of freedom within dense media, where heavy quark 
energy loss is carried out. 


studied in greater detail than in several earlier efforts. 
In this work, we have considered a hard virtual photon 
scattering off a hard heavy quark (within a proton), that 
converts it to a slow moving heavy quark that moves 
through the remainder of the nucleus before escaping and 
fragmenting into a jet containing a heavy meson. 

In this work both transverse broadening as well as the 
longitudinal drag and longitudinal diffusion, have been 
studied on an equal footing. We have categorically fo¬ 
cussed our study on “semi-hard” quarks where the mass 
and momentum scale as M,p ^ a/AQ, as these are the 
quarks for which mass modifications is most prominent. 
We have used power counting arguments loosely based on 
Soft Collinear Effective Theory (SCET) at various stage 
to isolate the leading contributions. It was shown in our 
earlier studies that both longitudinal and transverse mo¬ 
mentum transfers have a comparable effect on the off- 
shellness of the heavy-quark [32j |. This earlier work im¬ 
plied that longitudinal transfers, not only lead to the drag 
and diffusion, similar to light flavors, but will also notice¬ 
ably affect the radiative loss and left strong indications 
that for heavy quarks, the drag induced radiation may 
be as significant as transverse momentum diffusion (q) 
induced radiation. 

In this paper we have explicitly demonstrated that the 
gluon bremsstrahlung spectrum from a semi-hard heavy 
quark is indeed strongly modified by drag induced radia¬ 
tion. We have shown that due to the presence of the (—)- 
light-cone momentum exchange from the medium (fc“), 
in our calculations, the definition of all the transport co¬ 
efficients for heavy quark is different from that for light 
quark. Thus transport coefficients may indeed depend on 
properties of the probe i.e. mass or on the Whether 
this is phenomenologically significant cannot be ascer¬ 
tained at this point, and is left for a future investigation. 
This explicit dependence on the (—)-light cone momen¬ 
tum was absent in the limit of —>■ 0 , assumed in several 
prior calculations. 

The implications of the present study on the phe¬ 
nomenology of me is under way. It this work we have 
shown that the gluon bremsstrahlung spectrum of heavy 
quark (unlike light quark) is parametrically sensitive to e 
which quantifies the amount of drag the moving quark ex¬ 
periences. This result can be used to estimate the value of 
this sub-leading non-perturbative jet transport parame¬ 
ter (e) from heavy flavor data of HIC experiments. These 
extra additive contributions may lead to an eventual so¬ 
lution of the heavy quark puzzle. We leave these for a 
future effort. 
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VII. APPENDIX 



FIG. 2. 


There are three different cuts (central, left and right) in Pig.[2] and their contributions are, 
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In Fig. [3] there are only central cut, with the contribution. 
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FIG. 3. 



FIG. 4. 


In Fig. [4] there are two different cuts for induced gluon radition, central and left. 
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In Fig. [5] there are two different cuts for induced gluon radition, central and right. 
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FIG. 5. 
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FIG. 6. 


In Fig. [6] there are again all three possible cuts for induced gluon radition, 
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■^33 = exp 
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There are also three cuts in Fig. [7], and the contributions are, 
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The contributions from the three cuts in Fig. [5] are, 
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FIG. 9. 


Contributions from two cuts from Fig.[S] are as follows, 
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FIG. 10. 
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Contributions of Fig. [10] are, 
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FIG. 11. 


Fig. HI] has one possible cut, 
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Fig- [12] also has single possible cut. 
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FIG. 12. 
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